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We propose possible ways of explaining the net charge event-by-event fluctuations in Au-fAu 
collisions at the Relativistic Heavy Ion Collider within a quark recombination model. We discuss 
various methods of estimating the number of quarks at recombination and their implications for the 
predicted net charge fluctuations. We also discuss the possibility of diquark and quark-antiquark 
clustering above the deconfinement temperature. 
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Fluctuations of the net electric charge of all particles 
emitted into a specified rapidity window have been pro- 
posed as a possible signal for the formation of deconfined 
quark matter in relativistic heavy ion collisions P, 
The argument at the basis of this proposal is that charge 
fluctuations in a quark-gluon plasma are expected to be 
significantly smaller (by a factor 3 — 4) than in a hadronic 
gas. Because the net charge contained in a given volume 
is locally conserved and can only be changed by parti- 
cle diffusion, thermal fluctuations generated within the 
deconfined phase could survive hadronization and final 
state interactions. Quantitative estimates of the diffu- 
sion of net charge showed that the survival of these fluc- 
tuations from an early stage of the collision requires a 
moderately large rapidity window 0]. 

The most widely used measure for the entropy normal- 
ized net charge fluctuations is the D measure ^ : 

D = 4((Ag)2)/iVeh, (1) 

where ((AQ)^) denotes the event-by-event net charge 
fluctuation within a given rapidity window Ay, and Nch 
is the total number of charged particles emitted in this 
window. For a free plasma of quarks and gluons Z) w 1, 
while for a free pion gas D w 4. For the comparison with 
experimental data a number of corrections for acceptance 
and global charge conservation must be applied to the ex- 
pression for D 4] . The relation of the ZJ-measure to other 
measures of net charge fluctuations has been discussed by 
various authors lEIEH- 

Several experiments have measured net charge fluctu- 
ations in heavy ion collisions at the CERN Super-Proton 
Synchrotron (SPS) and at the Relativistic Heavy Ion Col- 
lider (RHIC) in Brookhaven H H, [H [ill . The results 
for D are generally somewhat smaller than 4, but much 
larger than the value predicted for a free quark-gluon 
gas. For example, the STAR collaboration has measured 
D = 2.8 ± 0.05 in central Au-f Au collisions at -^/sjvw = 
130 GeV 0, before applying corrections for global charge 
conservation and other effects The PHENIX exper- 
iment measured net charge fluctuations in a limited az- 
imuthal acceptance window around midrapidity, which 
extrapolate to a value Z) « 3 0. These results are sur- 
prising, because many other observables indicate that a 



deconfined quark-gluon plasma is formed in these colli- 
sions. 

Bialas has argued that the measured values of D could 
be compatible with the net charge fluctuations in a de- 
confined quark phase, if hadronization proceeds accord- 
ing to simple valence quark counting rules 12] and if glu- 
ons do not play an active role in the hadronization. In- 
deed, hadron abundances measured in relativistic heavy 
ion collisions at the SPS and RHIC are well described by 
combinatorial quark recombination models, such as AL- 
COR We here pursue this idea further and explore 
various scenarios of valence quark recombination in order 
to better understand how the puzzle posed by the mea- 
sured value of D can be resolved. We also discuss the 
constraints on such a resolution from the measured final- 
state entropy and the second law of thermodynamics. 

The recombination of thermalized valence quarks has 
recently been proposed as the dominant mechanism for 
the production of hadrons with transverse momenta of a 
few GeV/c in Au-|-Au collisions at RHIC [HEIEES 
0, ^1 . The RHIC data have provided compelling evi- 
dence for this hadronization mechanism |2Cll2l| . Valence 
quark recombination explains the enhancement of baryon 
emission, compared with meson emission, in the range of 
intermediate transverse momenta (roughly from 2 to 5 
GeV/c), and it naturally describes the observed hadron 
species dependence of the elliptic flow in the same mo- 
mentum region in terms of a universal elliptic flow curve 
for the constituent quarks [l5| . 

When one wants to describe hadronization by quark 
recombination not only at intermediate momenta, but 
over the entire hadron momentum range, except for very 
large momenta where parton fragmentation is thought to 
dominate, entropy becomes an important consideration. 
The naive application of recombination can easily entail 
a violation of the second law of thermodynamics, because 
the number of independent particles decreases when va- 
lence quarks recombine into hadrons. Greco et al. [l^ 
have argued that this problem may be circumvented by 
including the decay of hadronic resonances, such as the 
p-meson, into the calculation of the entropy balance. We 
next discuss the entropy problem in a more comprehen- 
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sive manner by considering the entropy content in re- 
alistic models of the hadronic phase (the resonance gas 
model) and the quark phase (lattice-QCD calculations). 

The equilibrium entropy per particle is only a fixed 
constant (3.60 for bosons, 4.20 for fermions) for free 
massless particles. For particle with mass, the entropy 
per particle is a function of the particle mass m and the 
temperature T. For m/T > 3 a good approximation is 
S/N = 3.50 + mjZT . The inclusion of mass is impor- 
tant, since a large fraction of the hadrons created at the 
moment of hadronization is quite heavy - the average 
hadron mass at chemical freezeout is about 800 MeV/c^ 
in the absence of medium induced modifications of the 
hadron masses. The average value of the entropy per 
hadron therefore significantly exceeds the canonical value 
{S/N)o w 4 in thermal equilibrium. Including only the 
pseudoscalar and vector ground state meson nonets and 
the octet and decuplet baryons, one finds {S/N) ~ 4.72 
at T = Tc = 170 MeV; including all known meson and 
baryon resonances, one obtains (S/N) = 5.15. 

An estimate of the final-state entropy per unit ra- 
pidity produced in central y^SNN = 130 GeV Au-|-Au 
collisions has recently been derived from experimen- 
tal data (hadron yields, spectra, and source radii) by 
Pal and Pratt: dS/dy = 4450 ± 400 H^. Using the 
measured charged multiplicity of dNch/dy = 526 ± 
2(stat) ± 36(syst) 0, this value can be converted into 
an estimate of the final entropy per particle of S/N « 
{dS/dy)/{l.5dNci,/dy) « 5.64 ± 0.6, which is slightly 
higher, but agrees within error with the equilibrium value 
for the full resonance gas at chemical freezeout. Indeed, 
numerical simulations of the reactions among hadrons in 
the hadron gas phase |2^ show that the particle num- 
ber increases only slightly (by about 5%) due to inter- 
actions after hadronization, compared with the number 
obtained by letting all hadrons formed at chemical freeze- 
out decouple and decay without reinteraction. A mod- 
est increase in the value of S/N during the hadron gas 
phase can easily be understood as the effect of volume 
expansion, since the entropy per particle for fixed par- 
ticle number increases logarithmically with the volume: 
S/N ^S,jN + ln{V/V,q). 

At the same time, the entropy content of the quark 
phase is strongly reduced due to interactions near Tc. 
Recently, th^CP-PACS collaboration ^ and the Biele- 
feld group j25| have calculated the pressure and energy 
density at finite temperature and zero chemical potential 
on the lattice. The Bielefeld simulation was done keeping 
the ratio niq/T fixed, while the one by the CP-PACS col- 
laboration was performed with the ratio of pseudo-scalar 
to vector meson masses, mps/my fixed. The CP-PACS 
collaboration used temporal lattice sizes Nt = 4, 6 and 
Nf = 2 quark flavors; the Bielefeld group used only 
Nt ~ 4, but explored a range of flavor multiplicities 
Nf = 2, 2-1-1, and 3. In order to extract physical quanti- 
ties from lattice QCD calculations, extrapolations (to the 



thermodynamic limit, continuum limit, etc.) are manda- 
tory. Because of a finer lattice spacing at a given tem- 
perature, the CP-PACS simulation for TVt = 6 is closer 
to the continuum limit and may be slightly better suited 
for the purpose of extracting the entropy density near T^. 
We use these results here. However, we should keep in 
mind that all lattice data for thermodynamic quantities 
are still obtained with unphysically large quark masses. 

We list the obtained values of e/T'^, P/T^, e/esB, 
P/PsB, and s/ssB, at mps/my — 0.65 and 0.80 below 
in Tables HI and Hll e, P, and s denote the energy den- 
sity, pressure, and entropy density, respectively, and esB, 
PsB, and ssB, are their values for the free gas of mass- 
less quarks and gluons on the lattice used in the simu- 
lation. The statistical errors are shown only for s/sgB- 
Obviously, the size of the systematic errors are still sub- 
stantial. However, at the same time, it is obvious that 
the entropy density of the quark-gluon plasma is con- 
siderably suppressed with respect to the corresponding 
Stefan-Boltzmann (SB) value near Tc. 



TABLE I: mps/mv = 0.65 case. All calculations are per- 
formed with {Nt — 6). Statistical errors are shown only for 

s/ssB. 

T/Tc e/T^ P/T^ g/esB P/Psb s/ssb 



0.92 8.03 0.88 

1.04 11.91 1.28 

1.30 13.67 2.75 

1.61 14.44 3.35 

2.00 12.93 3.96 



0.489 0.172 0.413 ±0.051 

0.725 0.250 0.612 ±0.043 

0.833 0.536 0.762 ±0.038 

0.879 0.653 0.826 ± 0.038 

0.787 0.772 0.784 ±0.038 



TABLE 11: Same as Table HI except mps/mv = 0.80. 

T/Tc e/T^ P/T^ e/esB P/Psb s/ssb 

0.80 3.73 0.17 0.227 0.033 0.181 ±0.054 

0.89 1.91 0.26 0.116 0.061 0.101 ±0.041 

1.12 12.08 1.62 0.736 0.316 0.636 ± 0.044 

1.38 11.98 2.66 0.730 0.519 0.679 ± 0.042 

1.67 11.80 3.54 0.719 0.690 0.712 ±0.039 



From these results we conclude the following: 

1. The entropy per particle in the hadronic gas, and 
therefore the entropy content of the hadronic phase 
at chemical freezeout, is considerably larger than 
often assumed. 

2. The entropy density of the quark phase is signifi- 
cantly suppressed near Tc, most likely due to cor- 
relations among the quasiparticles caused by their 
strong interactions. 
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These two conclusions make the recombination picture 
of hadronization more compatible with the entropy con- 
straint. Namely, if the quark-giuon plasma at hadroniza- 
tion consists of strongly interacting quasi-particles (pos- 
sibly constituent quarks) with strong correlations, and if 
many of the hadrons created at hadronization are heavy, 
quark recombination and the concomitant particle num- 
ber decrease can be more easily reconciled with the sec- 
ond law of thermodynamics. 

As we have seen, the increased entropy per hadron in 
the massive resonance gas may well allow for an isen- 
tropic transition from the deconfined to the confined 
phase by (sudden) quark recombination. At present, 
however, we cannot directly compare the entropy con- 
tent of both phases, since the volume at hadronization is 
not unambiguously known. Therefore, we cannot convert 
the entropy density determined on the lattice into a total 
entropy. Furthermore, lattice calculations do not tell us 
the number of (quasi-)particles, because there is no lat- 
tice definition of particle density. Therefore, a more de- 
tailed comparison of the entropy content of the hadronic 
phase and the quark phase at hadronization remains as 
a problem for future investigations. 

We now return to the calculation of net charge fluctua- 
tions in a bulk recombination scenario. The fluctuations 
of the net charge Q = qiUi are given by 



^cfk^n,){nk)q,qk (2) 



i,k 



(2) 

where c^^. are the normalized two-particle correlation 
functions: 
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(3) 



(4) 



The last expression in eq. Q shows that c^^."* is positive if 
there is a positive correlation between the quarks of fla- 
vors i and k. In the absence of two-particle correlations, 
(0) can be rewritten as: 

(<5g2) = liNu+ Nu) + ^{Nd + + N, + Ns) , (5) 

where Ni = (m) denotes the average number of con- 
stituent quarks of flavor i. 

Our strategy is now as follows: Knowing the number 
of final-state charged hadrons within a given rapidity 
interval, dNch/dy, we can extrapolate by means of the 
statistical hadronization model |2fiL to the thermal 
abundances of hadrons produced at the critical temper- 
ature Tr. We can then determine the total number and 



flavor distribution of valence quarks contained in these 
hadrons. Assuming valence quark recombination, using 
eq. (O, and neglecting correlations, we can then calcu- 
late the expected net charge fluctuation at hadronization. 
The prediction for {6Q'^) derived in this way can then be 
compared with the measured value of this quantity. 

We start from the measured charged multiplicity in 
central Au-|-Au collisions at 130A-GeV mentioned above. 
Including all established meson and baryon resonances 
in the chemical freezeout, we find that a total hadron 
number dN\^a_d/dy = 507 ±35 is needed at Tc to generate 
this charged particle number, when all sequential decays 
of unstable hadrons are taken into account. These 507 
hadrons contain a total of 1125 quarks and antiquarks. 
Applying eq. this gives a prediction 



d{SQ'^)q/dy = 286 ±23. 



(6) 



The measured value of the hadronic net charge fluctua- 
tion is lal 



d{6Q^)he,d/dy 



D X dN^Jdy 

X (2.8 ± 0.05) X (526±2±36) 



368 ± 33 . 



(7) 



We note that the errors in eqs. © and {Tj) are strongly 
correlated because they derive, in part, from the same 
uncertainty in the measured value of the charged particle 
multiplicity. Clearly, the net charge fluctuations resulting 
from the recombination of quarks derived in this way can 
explain only about 80% of the observed fluctuations. 

If, instead, we had included only the ground state 
nonet mesons and octet and decuplet baryons at chemi- 
cal freezeout, the inferred number of hadrons at Tc would 
have been larger, diVhad/c^y = 635 ±43, because fewer de- 
cays contribute to the final multiplicity. The number of 
recombining quarks would be accordingly larger, about 
1377, yielding the increased prediction 



d{6Q^)Jdy ^ 345 ± 29 , 



(8) 



which is close to the observed value eq. (0. 

A third way of estimating the number of constituent 
quarks at the moment of recombination is to start from 
the value dS/dy w 4450 derived by Pratt and Pal and 
use the calculated entropy per hadron [S/N — 5.15) for a 
resonance gas at to obtain an estimate for the number 
of hadrons at freezeout: dN\^^^/dy — 864±78. This value 
translates into an estimate for the net charge fluctuations 
from quark recombination of 



d((5Q2)q/dy = 487±39, 



(9) 



which lies significantly above the observed value. Con- 
sidering the systematic uncertainties inherent in these 
estimates, we may conclude that the observed net charge 
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fluctuations in Au+Au collisions at RHIC are compatible 
with the mechanism of bulk hadronization via recombi- 
nation of valence quarks in the absence of significant net 
charge correlations among the quarks. However, a more 
detailed estimate on the possible production of entropy 
in the hadronic phase would be desirable to better con- 
strain the analysis. Another source of uncertainty in our 
analysis is the possibility of the net-charge fluctuations 
increasing modestly in the hadronic phase due to diffu- 
sion 

We now discuss a modified variant of the recombina- 
tion process. Recently, quenched lattice QCD calcula- 
tions have shown evidence for the existence of mesonic 
bound state correlations even above the critical temper- 
ature 0, H H 0. Brown et al. argued within an 
effective field theory that bound states of charmed quark 
mesons, light quark mesons and gluons exist above 
[3^ . These findings suggest that qq and qq pairs may 
participate in hadronization mechanism as "elementary" 
constituents, just like individual quarks and antiquarks. 
In order to explore such a scenario, we modified eq. (0 
as follows: 



\ ~ qj] 



(10) 



where riij and the number of qq and qq pairs, 

respectively. For simplicity, we assume that the average 
number of qq (qq) pairs is proportional to the products of 
the individual quark numbers: (riij) = aNiNj] {uij) — 
f3NiNj, where a and (3 are the relative pairing weights. 
We have again neglected the correlation terms. 

The first term in eq. yields eq. The second 
term, which denotes the contribution from diquarks, is 
given by 



16 4 

= —aN^N^ + -a{NaNd + N,N, + NdNs) 

+ ^a{N^Nd + N^N,) + {q^q), 



(11) 



while the third term, denoting the contribution from 
quark-antiquark pairs, is: 

= P{NdNu+N,Nu + NuN^ + NuN,). (12) 

We can now constrain the parameters a and /3 using 
experimental value of (SQ^) for the two cases discussed 
above, counting only ground state hadrons or all known 
resonances at chemical freezeout. The total number of 
quarks and antiquarks on the right-hand side of eq. Ijlll) 
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FIG. 1: (a) Relation between the weight of the contributions 
from qq and qq pairs in the recombination process. The solid 
(dashed) line corresponds to a fixed total number of quarks 
and antiquarks of 1125 (1377). (b) Relation between the num- 
ber of individual quarks (antiquarks) and the number of di- 
quarks or qq pairs participating in the recombination process. 
Uq is the number of quarks and antiquarks contributing to the 
first term of eq. lllll . 



is constrained to be the same as that obtained from the 
statistical model. 

Figure 1(a) shows the relation between the weights of 
qq and qq pairs. This calculation is done in the simplest 
case {Nu = Na = Nd ^ N^, = Ng). The impor- 
tant point is the existence of the region where both a 
and /3 are positive, which confirms the possibility of a 
contribution from qq and qq pairs in the hadronization 
mechanism. We observe a similar tendency in both cases, 
counting ground state hadrons or all resonances. Diquark 
pairs are more favored than qq pairs in the sense that a 
solution with /? = is possible, but one with a = is 
not. In fact, the value of (3 is not well constrained by the 
charge fluctuations, because the difference between the 
charge of qq and qq in eqs. (|ll|l and (|12|) . For example, 
in the simplest case, J2ijili + ^ ^^-^u ^^'^ 

~ ^ 2/9 A^^, which implies that qq pairs 

are favored by a factor of about 2. Perturbative QCD 
suggests that the qq channel is more attractive than the 
qq channel, implying a < (3 m eq. (|ll|l . In Fig. 1(a) the 
lower part of the dashed line and the whole solid line are 
forbidden if this condition is imposed. However, since 
the lattice results indicate that hadronization occurs via 
strong interactions between plasma quasi-particles, it is 
not clear that this perturbative argument is applicable. 

In order to clarify the relative numbers of diquarks, 
quark-antiquark pairs, and individual (anti-)quarks par- 
ticipating in the recombination process, we plot the re- 
lation among A'^ for quarks and antiquarks, Nij 
for diquarks, and Nij for quark-antiquark pairs in 
Fig. 1(b). In the full resonance gas scenario, for example, 
we have Y.^ = 134, J2^J % = 493, and J2^J N^J = 2, 
showing that qq clustering dominates. The numbers 
of quarks and antiquarks decrease linearly as those of 
diquarks and quark-antiquark pairs increase along the 
solid and dashed lines, respectively. In the ground state 
hadrons only scenario, we find a region where recombi- 
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nation from individual quarks and antiquarks is domi- 
nant. On the other hand, in the full resonance gas sce- 
nario, hadrons are predominantly created from diquark 
or quark-antiquark clusters, which is difficult to recon- 
cile with the elliptic flow data from RHIC j^^, 0| , since 
that data strongly suggests a constituent quark counting 
rule El E [IE El- If we demand that the weight of 
individual quarks and antiquarks in the recombination 
process should be so large that the quark number scaling 
approximately survives, the full resonance gas scenario is 
disfavored. 

In summary, we have investigated charged particle fluc- 
tuations at RHIC in the framework of the parton recom- 
bination model of hadronization and find that within the 
present systematic uncertainties parton recombination is 
compatible with the measured charged particle fluctua- 
tions. We found that the behavior of the entropy density 
for an interacting deconfined system close to Tc and the 
entropy per particle for the massive resonance gas sup- 
port the recombination picture. Finally, we have investi- 
gated the possibility of bound state correlations above Tc 
and find them consistent with the parton recombination 
approach as well, albeit constrained by the valence quark 
number scaling observed in the data. 
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